FREE DETERMINISTIC EQUIVALENTS, RECTANGULAR 
RANDOM MATRIX MODELS, AND OPERATOR- VALUED FREE 

PROBABILITY THEORY 



ROLAND SPEICHER AND CARLOS VARGAS 
WITH AN APPENDIX BY TOBIAS MAI 



Abstract. Motivated by the asymptotic collective behavior of random and 
deterministic matrices, we propose an approximation (called "free determin- 
istic equivalent") to quite general random matrix models, by replacing the 
matrices with operators satisfying certain freeness relations. We comment on 
the relation between our free deterministic equivalent and deterministic equiv- 
alents considered in the engineering literature. We do not only consider the 
case of square matrices, but also show how rectangular matrices can be treated. 
Furthermore, we emphasize how operator-valued free probability techniques 
can be used to solve our free deterministic equivalents. 

As an illustration of our methods we show how the free deterministic equiv- 
alent of a random matrix model from [6] can be treated and we thus recover 
in a conceptual way the results from [5]. 

On a technical level, we generalize a result from scalar valued free probabil- 
ity, by showing that randomly rotated deterministic matrices of different sizes 
are asymptotically free from deterministic rectangular matrices, with amalga- 
mation over a certain algebra of projections. 

In the Appendix, we show how estimates for differences between Cauchy 
transforms can be extended from a neighborhood of infinity to a region close 
to the real axis. This is of some relevance if one wants to compare the original 
random matrix problem with its free deterministic equivalent. 



1. Introduction 

Voiculescu's results [22l [23l [TTl [15] on the asymptotic freeness of random ma- 
trices have provided a way to understand the collective asymptotic behavior of a 
variety of random matrix ensembles. Many known results from the random matrix 
literature have found straight-forward solutions via free probability. Additionally, 
this perspective has lead to new results (e.g. the asymptotic freeness of randomly 
rotated deterministic ensembles). The engineering community, in particular, has 
added concepts and results from free probability theory to their toolbox for deal- 
ing with random matrix ensembles in the context of wireless communications, see, 
e.g., [5D|. However, still there is the widespread perception that free probability 
cannot deal with situations where rectangular matrices are involved or where one 
does not have an asymptotic eigenvalue distribution of some of the involved ensem- 
bles. In the latter context, so-called "deterministic equivalents" (apparently going 
back to Girko [5], see also [3]) have recently been studied quite extensively in the 
engineering literature. 

In this paper we will show 

• how those deterministic equivalents can be understood and formalized in a 
conceptual way by using free probability theory 

• and that this can actually also be done in the presence of rectangular ran- 
dom matrices 
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Operator- valued free probability [HI (TH] will play a crucial role in these investi- 
gations. Firstly, it is actually needed to deal in a nice way with rectangular random 
matrices - following the ideas of Benaych-Georges @] we will embed the rectangular 
matrices in a large matrix space, which has the structure of a projection- valued 
probability space (which is one of the simplest cases of operator- valued spaces). 
Secondly, for more advanced random matrix models one usually has to rely on 
operator-valued freeness results to derive the equations for the deterministic equiv- 
alents. The fundamental observation that the operator valued version of free prob- 
ability allows to understand a broader class of matrix ensembles, such as those 
consisting of band matrices or block matrices, goes back to Shlyakhtenko [TB]. See 
also [T7] for a first use of operator-valued freeness in the context of wireless com- 
munication problems. 

In order to illuminate, and also to show the power of our general approach we will 
re-derive the equations for the deterministic equivalent of the Cauchy transform of 
a model considered in [5] (in the asymptotic regime this model was also considered 
before in [H]). The model is the following: 

(1.1) $^ = )*, 

i=l 

where is a fixed integer and 

(Al) For each I <i < k, {N^^'')Ni^fq is a sequence of natural numbers such that 

N^^^/N^c, e (0,oo). 
(A2) For each 1 < i < fc, H^'^'^ is an x N^^'' deterministic matrix and Tj*-^-* is 

an N^^^ X N^^^ deterministic matrix. 
(A3) The matrices U^'^'^ are independent Haar-distributed unitaries bomU{Nj^^^). 
The paper is organized as follows: 

Section 2 is devoted to the basic theory on operator-valued free probability. 
We state some results from [1], generalizing the asymptotic freeness results by 
Voiculescu to the case where we allow rectangular matrices of different sizes. These 
results are stated in the framework of rectangular probability spaces. Our con- 
tribution to the theory is a generalization of the asymptotic freeness of randomly 
rotated deterministic matrices. 

In Section 3 we define the notion of a "free deterministic equivalent" for a ran- 
dom matrix model, based on our understanding of the limiting collective behavior 
of random and deterministic rectangular matrices provided in the previous section. 
We will also point out the relation between our free deterministic equivalent and 
the deterministic equivalents from the engineering literature: whereas the latter 
is an approximation on the level of equations for the Cauchy transforms, our ap- 
proximation is directly on the level of operators; but in the end they coincide, 
i.e. the Cauchy transform of our free deterministic equivalent satisfies exactly the 
deterministic equivalent equations. 

In Section 4 we present the free deterministic equivalent of the random ma- 
trix model (jl.ip and derive, by using operator-valued free probability results, the 
equations for its Cauchy transform. This will recover the results from [5] . 

In the last section we give the proof of our result on randomly rotated rectangular 
matrices from Section 2. 

In the Appendix (written by T. Mai) we show how estimates for differences 
between Cauchy transforms can be extended from a neighborhood of infinity to a 
region close to the real axis. This is of some relevance if one wants to compare the 
original problem with its free deterministic equivalent. 
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2. Operator- Valued Free Probability 

2.1. Preliminaries. We begin by recalling the basic concepts on operator-valued 
free probability, for more details see [211 [T9] . 

Let B he a unital algebra. A B-probability space is a pair {A, F) consisting of an 
algebra A containing B and a linear map F : A B satisfying 

F (bab') = bF{a)b', V6, b' eB,aeA 

and 

F(l) = l. 

Such a map F is called a conditional expectation. From the properties above we 
observe that F (&) = 6 for all & G ;B. 

A collection of subalgebras B d Ai, . . . , An C ^ is free with amalgamation over 
B (or, B-free, for short) if we have that for every /c G N 

F (ai(i) . ..a^k)) = 

whenever we have F (oi(j)) = and aj^) G ^^(j) for some e{j) G {!,..., n}, 
1 < j < and that e (1) ^ e (2) =^ ■ • • / e (k). (Note that e (1) = e (3) is allowed.) 

Elements oi, . . . , a„ are B-free if the algebras {ai,B) , . . . , (a„, B) are S-free. 
{{ai,B) is here the subalgebra of A which is generated by oi and by B.) 

For n G N, a C- multilinear map / : A"" B is called B -balanced if it satisfies 
the S-bilinearity conditions, that for all b,b' G ai,...,a„ G A, and for all 
r = 1, . . . , n — 1 

/ (5ai, . . . ,a„5') = 6/ (ai, . . . , a„) 5' 

/ (ai, . . . , a^b, a^+i, . . . , a„) = / (oi, . . . , a^, foa^+i . . . , a„) 

For a finite totally ordered set iS, a partition tt = {Vi, . . . , T4} is a decomposition 
of S into pairwise disjoint non-empty subsets, called blocks. A partition is non- 
crossing if whenever we have a<b<c<d£S such that a,c ^ Vi and b,d € 
Vj, then Vi = V, . We denote by NC{n) the set of non-crossing partitions of 
[n] := {1, . . . , n} and NC :— lJn>o NC{n). The partial order of reverse refinement 
(cr < TT iff every block of a is completely contained in a block of tt) turns NC{n) 
into a lattice. We denote the smallest and largest partition of NC{n) by 0„ and 
1„, respectively. 

A collection of ;B-balanced maps (/7r)^g7vc ^^^"^ multiplicative (w. r. t. 

the lattice of non-crossing partitions) if, for every tt G A^C, is computed using 
the block structure of tt in the following way: 

• If TT = 1„ G NC (n), we just write /„ :— /^. 

• If 1„ 7^ TT = {Vi, . . . , Vfe} G NC {n) , then by a known characterization 
of NC, there exists a block Vr = {s + 1, . . . , s + 1} containing consecutive 
elements. For any such a block we must have 

/tt (ai, ■ ■ ■ , a,n) = /7r\K ('Jli ■ ■ ■ : 0,sfl (Os+l, ■ ■ ■ , O-s+l) , ^g+l+l, . . . , a„) , 

where 7r\14. G NC {n — I) is the partition obtained from removing the block 

Vr. 

We observe that a multiplicative family {fn)^^iqc entirely determined by 
ifn)n(Efq- On the other hand, every collection (/n)„gpj of S-balanced maps can 
be extended uniquely to a multiplicative family (/7r)^gjvc- 

The operator-valued cumulants {k,-k)t^^i<ic ^'^^ indirectly and inductively defined 
as the unique multiplicative family of Z?-balanced maps satisfying the (operator- 
valued) moment-cumulant formulas 

F(ai...a„)= ^ (ai, . . . ,a„) 

TTeNC(n) 
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The conditional expectation yields also a multiplicative family by simply defining 
F„ : ^" B by 

F„ (ai, . . . ,a„) = F (ai . . . a„) , 

and extending multiplicatively to NC. Actually, the relation between cumulants 
and moments is much richer. The cumulants can also be obtained from the moments 
via a Mobius inversion in the lattice of non-crossing partitions, and one has, for 
example 

(ai, . . . ,a„) = ^ Fcr (oi, . . . ,a„)^„[tT,7r], 

(T<7r 

where /i„ : NC {nf' — > C is the Mobius function in NC (n). 

By the cumulants of a tuple ai, . . . , afe ^ A, we mean the collection of all cumu- 
lant maps 

(61, . . . ,fe„_i) H> (0^161,0^262, ••• ,ai„) 

for rt G N, 1 < < k. Wc sometimes write ^f^'^i'v'"'' to emphasize the 

underlying subalgebra. 

A cumulant map is mixed if there exists r < n such that v 7^ V+i- 

The main feature of the operator-valued cumulants is that they characterize 

freeness with amalgamation [19| : The random variables oi, . . . , a„ are ;B-free iff all 

their mixed cumulants vanish. 

The operator-valued Cauchy transform is given by 

^ ^ ri>0 

The operator-valued TZ-transform is given by 

^?W = E^" (0,6a,..., 6a). 

n>l 

The vanishing of mixed moments for free variables implies the additivity of the 
cumulants, and thus also the additivity of the 7?.-transforms [21j : If oi and 02 are 
S-free then we have for 6 G S that TZai+a2 (b) ~ Ti-ai (6) -I- TZa^ (6). 

As in the scalar case, these transforms satisfy the functional equation 

Gf(6) = (7^« (Gf(6))-6)-^ 

The above transforms and the functional equation are on one hand just identities 
for formal power series, but one can also consider them as relations for analytic 
functions, on appropriate domains. For more on the analytic properties in the 
operator- valued context, see pTl [2]. 

2.2. Working on Different Levels. One of the most fascinating aspects about 
the operator- valued version of free probability is that one can consider several condi- 
tional expectations on a single space, and compare the corresponding distributions. 

In order to guarantee the existence of nice conditional expectations onto subal- 
gebras one needs some analytic structure. Particular nice is the case of a tracial 
W*-probability space {A, r), which means that ^ is a W*-algebra, or von Neumann 
algebra, (i.e., a unital subalgebra of bounded operators on a Hilbert space, closed in 
the weak operator topology), and that r is a normal state which satisfies the trace 
condition r(aia2) = T(a2ai) for all 01,02 G A. We will only encounter situations 
of this type. In particular, matrices (equipped with the usual trace as state) and 
their limits are of this type. 
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If we start from a tracial H^*-probability space {A, r), then every VF*-subalgebra 
C C B C A induces a unique conditional expectation F : A ^ B compatible with 
r in the sense that t — t oF. 

This means in particular that the scalar-valued Cauchy transform can be ob- 
tained by the projection of the operator-valued one, namely 

When dealing with elements ai, . . . ,a„ in A which are not free, very often we 
can "enlarge" the algebra of scalars C to a subalgebra C C ;B C ^, in such a 
way that ai, . . . , a„ become S-free. Then, a variety of results can be used to treat 
and simplify the ^-distribution of products and sums of these free elements. The 
scalar-valued distribution can then be obtained in the end by projecting with r. 

We recall two results from [M] . The first proposition gives conditions for (operator- 
valued) cumulants to be restrictions of cumulants with respect to a larger algebra. 

Proposition 2.1. Let lEVcBcAbe algebras such that {A, E) and {B, F) are 

respectively B-valued and V-valued probability spaces (and therefore {A, F o E) is a 
V-valued probability space) and let ai, . . . , G A. 

Assume that the B-cumulants of ai, . . . , ak & A satisfy 

for all n Cz 1 < ii, . . . , i„ < fc, di, . . . , c?„_i G V. 

Then the V-cumulants of ai, . . . , are just the restrictions of the B-cumulants 
o/ ai, . . . , Ofc, namely 

«fi'!.'.;i„'''' {di,---, dn-i) = ^f^'^'Z''"'' {di, . . . , dn-i) , 
for all n Q N, 1 < ii, . . . ,in < k, di, . . . , c?„_i G "D. 

The second proposition gives a characterization of operator-valued freeness by 
the agreement of operator-valued cumulants in different levels. 

Proposition 2.2. Let 1 e V C B,J\f C A be algebras such that {A,F) and {B,F) 
are respectively B-valued and V-valued probability spaces. Then the first of the 
following statements implies the second. 

i) The subalgebras T>,M are free with amalgamation over B. 

ii) For every fc G N, ni, . . . , rife G J\f, 6i, . . . , bk-i G B, we have 

. . . ,?ifc_i6fe_i,nfe) = {niF , . . . ,7ife_iF (bk-i) ,nk.) . 

Moreover, the two statements become equivalent if we add the faithfulness condition 
onF : B ^V, that i/F (6162) = for all 62 e B, then hi = 0. 

2.3. Rectangular Spaces. Let [A, r) be a tracial Vl^*-probability space endowed 
with pairwise orthogonal, non-trivial projections pi, . . . ,pk adding up to one. Let 
V := (pi, . . . ,pk) denote the W^*-algebra generated by {pi, . . . ,Pk}- Then there 
exists a unique conditional expectation F : A ^ T> such that t o F = r. We 
actually have that 

fc 

(2.1) F(a) = Y,P^r{Pi)~'T{p,a). 

i=l 

With this, {A, F) becomes a 2?-valued probability space. 

These kind of projection-valued spaces are called rectangular probability spaces. 
We will denote by A'^^'^^ the set of elements a € A such that a = Piopj. Elements 
in Ui<i j<fc called simple and we write A^'^^ := A^^'^K Very often we will 
be interested in the compressed spaces {A'^^\t'^^^), where t^*^ (a) = T{pi)~^T (a) , 
for a G A^'\ 
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These spaces were introduced by Benay eh- Georges in [3J|3]. They allow to treat 
rectangular matrices of different sizes in a single space. 

From now on, we will restrict the use of the caligraphic letter T) to the W*- 
algebra generated by the projections of a rectangular probability space. Therefore, 
I?-freeness will always mean freeness with amalgamation over the algebra of pro- 
jections which defines the given rectangular probability space. It will be useful to 
specialize the corresponding notion of asymptotic freeness as well. 

Definition 2.3. Let {A, r) he a {pi, . . . ,pk) -rectangular space and let (.A„,r„)„>i 
be a sequence of (pi""* , ■ . ■ , p^"^ ) -rectangular spaces. Let oi , . . . , G A and a^"\ . . . , Um^ € 
An be collections of simple elements. We say that (a^"\ . . . ,am^) converges in T>- 
distribution to (ai, . . . , Om) • ■ • , o,m\pi^\ ■ ■ ■ jpI"'') converges in distribution 

to (ai, . . . , am, Pi, ■ ■ ■ ,Pk), arid we write 

\a\ , . . . , a\^' I (fli, . . . , Om) as n ^ oo. 

// oi, . . . , Om are (pi, . . . ,pk)-free, we say that a^"\ . . . , al^J are asymptotically V- 
free. 

Example 2.4. Consider the matrices from (II. ip for a fixed A^. We set Nq := N and 
let 

i=0 

We consider pairwise orthogonal projections Po,...Pk in the space {Am,tm) of 
M X M random matrices, where each Pi is a projection to a subspace of size 
Ni. Then we can embed our matrices Ui,Ti, Hi, Hi . . . ,Uk,Tk, Hk,H^ as simple 
elements in Am, as illustrated below: 



Po 


Hi 




Hk 




Ti,Pi 




H* 






Ui,U*i 














Tk,Pk 


m 








Uk,Ul 



We denote by A the embedding of A. A simple element in a (random) matrix 
rectangular space will be called a simple matrix. 

In order to achieve asymptotic freeness we need some requirements on the as- 
ymptotic distribution of the involved deterministic matrices (which we try to keep 
to a minimum). We will assume: 

(A4) Each matrix Ti converges in distribution w.r.t. ]^Tr, and HiH^, . . . , HkH^ 
converge in joint distribution w.r.t. ;^Tr. 

On these rectangular spaces, we are prevented from multiplying matrices that do 
not fit together (e.g. HiHj = 0). These kind of mixed moments do not show up 
in the distribution of so we will still be able to calculate this distribution by 

working in these rectangular spaces. The advantage here is that even when we are 
only asking for the existence of the limiting joint moments of HiH*, . . . , HkH'^, we 
will actually have that Hi , . . . , Hk have a limiting joint ^-distribution. 
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We recall results by Benaych-Georges in [3] generalizing, in the framework of 
rectangular spaces, Voiculescu's results on asymptotic freeness of square matrices. 

Proposition 2.5. Let k > 1 be fixed and for each N > 1, let {An,tn) be a 
(P^^\ . . . , pjf^^ ) -rectangular probability space of random matrices, such that tn(P^^'^) 
Ci G (0, 1], ? = 1, . . . , fc. Let ([/j'^'*)i>i be a collection of independent simple random 
matrices in An, such that for every is a Haar- distributed unitary 

ensemble in the compressed space (-Aj^j*'*'''', t^^*''"') for some I < j (i) < k. 

Furthermore, let Dn = (-Dj-^'')i>i be a collection of simple deterministic ma- 
trices, such that for all N, D^^^ G for some 1 < r{i),s{i),< k, and 
assume that (-Dj-^'')i>i converge in D- distribution. 

Then Dn, u[^\ U2^\ . . . are asymptotically D-free. 

Remark 2.6. As should be expected, if in the preceding theorem one considers 
Gaussian or non self-adjoint Gaussian simple matrices instead of Haar unitary sim- 
ple matrices, the same P-freeness results hold, where the corresponding random 
matrices converge to semicircular or circular simple elements in a limiting rectan- 
gular space. 

In [3] the case where rjv [Pi) -t- is also addressed. 

From the combinatorial point of view, the computation of the asymptotic scalar- 
valued moments of a product of rectangular matrices is performed in the same way 
as in the usual, square case, namely, it can be obtained as a sum of products of 
moments of smaller degree, running over certain collections of non-crossing parti- 
tions. These phenomena will be better observed later in our proof of the asymptotic 
freeness of randomly rotated deterministic matrices. 

Roughly speaking, the only difference from the square case is that one needs to 
scale each moment by taking the different sizes of the matrices into consideration. 
This scaling is automatically performed by the use of the ^-distribution, which 
puts the right weights on each moment. In particular, the order of convergence to 
the moments of free operators is the same as in the square case. 

We also want to remark that, by following the arguments for proving asymptotic 
freeness between deterministic and Wigner matrices (see [I] [12]), the previous re- 
sults should admit a straight-forward generalization to the case where the Gaussian 
matrices are replaced by Wigner matrices. 

2.4. Randomly Rotated Matrices. A well-known result from scalar-valued free 
probability states that one can consider a pair of collections pj^-* , XJj^-* of determin- 
istic matrix ensembles, where each collection of ensembles converges in distribution 
(with respect to -^Tr). Then if we randomly rotate all the matrices in one of the 
collections by conjugating with a Haar-distributed unitary matrix, the resulting 
collections 'D[^\UnT^2^^U^ are asymptotically free. 

We do not require the existence of an asymptotic joint distribution of v[^^ , . 
It follows from the Weingarten formulas for the joint distribution of the entries of 
Haar-distributed matrices, that these mixed moments do not appear in the calcu- 
lations, see, e.g., [51 [TS]. 

The following theorem is a generalization of the previous result in the rectangular 
framework. 

Theorem 2.7. Let k > 1 be fixed and for each N G N, let (^Ar,Eo ;^Tr) be a 
space of NxN random matrices with the structure of a {P^'^\ . . . , P^'^^) -rectangular 
probability space, with simple elements U^^^ G A^j^ such that: 
i) limAr^oo^Tr(/^^'^)) =00, 
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ii) Ui^\ . . . , ujf^^ are independent (entrywise). Each Ui is a Haar- distributed 
unitary random matrix in the compressed space {A^j!^ , r^'"* ) . 

Put Ujv := U[^'' + ■■■ + Uf^ and let 

Vi""^ = . . . , C^f)} and ) = [d[''\ )} 

be collections of simple deterministic matrices, each one with asymptotic "D- distribution. 
Consider 

UNVi''^V*N ■■= {VnDV% : D e 2?f . 

Then V^^"^ and Ujv^?2^''U^ are asymptotically V-free. 

If we assume that 'd[^^ , 2?2^'' have a joint distribution, then the previous result 
is an easy corohary from Proposition 12.51 Without such an assumption, the proof 
of Theorem 12.71 fohows the hues of the square case, and we wih leave it to the 
last section. First we want to show how it enables us to compute the asymptotic 
distribution of from Equation (jl.lD . 

Corollary 2.8. Let f/}^', . . . , t/f \ h[''\ . . . , i/f ^ and T*^^ . . . , Tf ^ be deter- 
ministic ensembles satisfying conditions (Al)-(A^) and consider their embedings in 
the rectangular spaces described in Example \2.4\ Then 

. . . , \ )c/f )*, . . . , c/rri^)?7f 

converges in joint distribution to a collection of elements (/ii, . . . ,hk,ti, . . . , tu) in 
a W* -probability space {A, t), which has the structure of a {po, . . . ,pk) -rectangular 
probability space, satisfying 

i) T(pi) = Ci. 

ii) and ti 

iii) {hi, . . . , hk) ,ti, . . . ,tk are free with amalgamation over {po, . . . ,Pk) ■ 

We could use the previous corollary to describe the asymptotic distribution of 
$Ar. Instead, we will introduce in the next section an approach that allows us to 
work with finite N , motivated by our knowledge of the limiting behavior of the 
involved matrices. 



3. Free Deterministic Equivalents 

To simplify our exposition, we first introduce the free deterministic equivalents 
for the case where all the matrices are square and have the same size. Then it will 
be natural to generalize it to the case of rectangular matrices with different sizes. 

3.1. Square matrices. Voiculescu's asymptotic freeness results on square random 
matrices state that if we consider tuples of independent random matrix ensembles, 
such as Gaussian, Wigner or Haar unitaries, their collective behavior in the large N 
limit is that of a corresponding collection of free (semi-) circular and Haar unitary 
operators. Moreover, if we consider these random ensembles along with determin- 
istic ensembles, having a given asymptotic distribution (with respect to the nor- 
malized trace) , the corresponding limiting operators become free from the random 
elements. 
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matrices: 








r yW 




Yiv = 




••'-^^(2) 


Uat = 




••'^j(3) 


Djv = 


(<),. 


• • ' ^i(4) 



We consider now a collection of independent random and deterministic N x N 

I : independent self-adjoint Gaussian matrices 

I : independent non self-adjoint Gaussian matrices 

I : independent Haar-distributed unitary matrices 

I : deterministic matrices, 

and a polynomial in non-commuting variables (and their adjoints) evaluated in our 

matrices 

Relying on the above mentioned asymptotic freeness results, we can then com- 
pute the asymptotic distribution of Pjv with respect to the expected normalized 
trace E o ;^Tr ~: t^, (which, for the case of self-adjoint polynomials, = Pn, 
coincides with the expected spectral distribution of Pjv) by going over the limit. 
We know that we can find collections S, C, U, D of operators in a non-commutative 
probability space {A,t), 



s 


= {si, ■ 




free semi-circular elements 


c 


= {ci,. 


••,C,(2)} 


free circular elements 


u 


= {ui,. 


• , '"i(3) } 


free Haar unitaries 


D 


= {di,. 


• • ,C?i(4)} 


abstract elements, 



such that S, C, U, D are *-free and the joint distribution of di, . . . (ii(4) is given by 
the asymptotic joint distribution of d[^^ , . . . , -D^^j . Then, the expected asymptotic 
distribution of Pn is that of 

P (si,. . . ,Si(i),Ci,. . . ,Cj(2), Wl, . . . ,Uj(3),di, . . . ,rfi(4)) =: Poo, 

in the sense that, for all k, 

Jim r;v(P^)=r(P^). 

In this way we can reduce the problem of the asymptotic distribution of Pn to 
the study of the distribution of Pqo , and this can be done by using the tools of free 
probability, due to the nice freeness relations exhibited by S, C, U, D. 

A common obstacle of this procedure is that our deterministic matrices may not 
have an asymptotic joint distribution. But now that we understand the limiting 
behavior of these random and deterministic matrices, it is natural to consider, for 
a fixed N, the corresponding "free model" 

P (^si, Si(i), ci, . . . , Cj(2), Ul, . . . , Uj(3), . . . , =: Pn 

where, just as before, the random matrices are replaced by the corresponding free 
operators in some space {An, Vn), but now we let the distribution of d[^\. . . , rfj^j 

be exactly the same as the one of d[^\ . . . , D^^J^ w.r.t ;^Tr. The free model P^ 
will be called the free deterministic equivalent for Pn- 

Once again we will be able to deal with P^ by the means of free probability. 
The difference between the distribution of P^ and the expected distribution of 
Pjv is given by the deviation from freeness of Xjv,Yjv,Ujv,Djv, the deviation of 
Xjv, Y;v from being free (semi)-circular systems, and the deviation of Ujv from a 
free system of Haar unitaries (for the case of Ujv, the matrices are already Haar 
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unitaries individually, but they are not free) . Of course for large N these deviations 
get smaller and thus the distribution of becomes a better approximation for 
the expected distribution of Pn. (By the concentration of measure phenomena, 
the difference between the empirical distribution of Pn and its expectation will 
typically also tend to zero in the large limit.) 

Let us denote by Gat the Cauchy transform of Pn and by the Cauchy 
transform of the free deterministic equivalent P^ . Then the usual asymptotic free- 
ness estimates show that moments of Pn are, for large A^, close to corresponding 
moments of P^ (where the estimates involve also the operator norms of the deter- 
ministic matrices). This means that for iV — >■ oo the difference between the Cauchy 
transforms Gn and G^ goes to zero, even if there does not exist individual limits for 
both Cauchy transforms. (In the appendix we will show how to extend estimates for 
the difference of Cauchy transforms in a neighborhood of infinity - corresponding 
to knowledge about the moments of the involved distributions - to a region close 
to the real axis; which can then be used to compare the respective distribution 
functions, or to estimate the Kolmogorov distance between the distributions.) 

In the engineering literature the notion of a deterministic equivalent (apparently 
going back to Girko 0, see also [9]) has recently gained quite some interest. This 
deterministic equivalent consists in replacing the Cauchy transform Gat of the con- 
sidered random matrix model (for which no analytic solution exists) by a function 
Gn which is defined as the solution of a specified system of equations. The specific 
form of those equations is determined in an ad hoc way, depending on the consid- 
ered problem, by making approximations for the equations of Gat, such that one 
gets a closed system of equations. 

In all examples of deterministic equivalents we know of it turns out that actually 
the Cauchy transform of our free deterministic equivalent is the solution to the 
equations of the deterministic equivalents, i.e., that Gn = G^. We think that our 
definition of a deterministic equivalent gives a more conceptual approach and shows 
clearly how this relates with free probability theory. In some sense this indicates 
that the only meaningful way to get a closed system of equations when dealing 
with random matrices is to replace the random matrices by free variables. We want 
to point out that the same phenomena was essentially also observed in |13j in the 
context of the so-called CPA approximation (a kind of mean-field approximation) 
for the Anderson model from statistical physics. In our present language their result 
can be rephrased as saying that the free deterministic equivalent of the Anderson 
model is given by the CPA approximation. 

Since in engineering applications rectangular matrices are quite typical, it is 
important to be able to have the notion of a free deterministic equivalent also for 
problems involving rectangular matrices. In the next section we will show how 
the above can be generalized to the rectangular setting, thus hopefully refuting 
the common perception that rectangular matrices cannot be dealt with in free 
probability. 

3.2. Rectangular matrices. In order to include models with rectangular matrices 
of different sizes, we have to consider slightly more complicated situations which 
take place in rectangular probability spaces. The main idea is that the space is split 
into pieces {A^j!f''^)i<ij<k and one can think of a rectangular matrix as an element 
in one of these pieces; in particular, square matrices will be regarded as elements 
in the compressed spaces A'"^ . 

We can still use all the results on asymptotic freeness for square Gaussian, Haar 
unitary and deterministic matrices in this framework, by just working in the corre- 
sponding compressed space {A^j^^ , r^'* ) , with the novelty that different sized square 
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random matrices are I?-free from each other and from deterministic rectangular 
matrices. 

To obtain the free deterministic equivalent, we will replace the random square 
matrices by free operators ((semi-)circulars and Haar unitaries) in the correspond- 
ing compressed spaces, just in the way described in the square case. When it comes 
to rectangular deterministic matrices in , we will replace them by operators, 
keeping the same distribution, but prescribing them to be D-free from the (semi-) 
circulars and the Haar unitaries. 

It is important to impose some restrictions in our non-commutative polynomials 
for this rectangular situation. We should only treat polynomials where the sums 
and the products are performed in such a way that the sizes fit. The reason for this 
is that our rectangular spaces prescribe non-compatible products of matrices to be 
zero, and hence, such products should not appear in our polynomials. 

Remark 3.1. In order to simplify notation, when the considered non-commutative 
polynomial P depends on a randomly rotated deterministic matrix [/(^)d(^^?7(^^* 
and not actually on C/^^^ D^^^ G A^^^ individually, we will just write (instead 
of udWu*) in our free deterministic equivalent, keeping in mind that d(^) win be 
free from all other deterministic matrices. 



4. An Example 

From our discussion on the previous section, the free deterministic equivalent for 
$ AT from (jl.ip will be 

i=l 

where (after supressing the N super-index) each hi and ti is a simple element in 
a (po; • ■ • 7Pfc)-pi'obability space {A,t), with conditional expectation F : A ^ V, 
satisfying 

i) ^(k) ^ffiorl<i<k 

ii) hi G A^°''\ U e A^'\ with 

(im)^i_Tr(i;'"), 
for 1 < i < fc, and for all to G N and all 1 < ii, . . . , im < fc we have 
t(") hi • • • h*^ ) = ^Tr {H,, HI--- H,^ H*^ ) 
t('i) {hi - - ■ h,^h*^h,,) = -^Tr [HI - - - H,^H*^H,,) 

iii) {hi, . . . ,hk} ,ti, . . . ,tk are P-free. 
We consider the subalgebras 

C := {V,hihl,...,hkhl) cA and C^°'> := poCpo C A^"'> , 

with the conditional expectations 

E:A->C and E^o' : ^ C^^^ 

compatible with t and t^^\ respectively. 

We are ultimately interested in the distribution of $ = $^ considered as an 
element of A^^K It is easy to see that E|_^(o) is C*-"' -valued and defines a conditional 
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expectation compatible with r^'^^ . Hence, by uniqueness, E^*^) must be E|_4(o) . Then 
we observe that for invertible c G C 

n>0 

= ^ E (poc^^Po (^Poc" Vo)") 

n>0 

= {ipocpor' (<& (Pocpor'Y) 

= GJ' (pocpo) , 

which follows from the commutativity of C with the elements of P (note that po € 
T>). Hence it suffices to describe the C- valued distribution of $. 

For this we first notice that B := {V, hi, ... , hk) and M := (P, ti, . . . , tk) are 
■D-free, so we are in position of applying Proposition [521 and we obtain 

= hi-^ K„ (ti-^ /ij^ bihi^ , ft-jj ^2^13 I ■ • ■ J tin ) 

= hi^K^ {U,F {h*^bihi^) ,ti^F {h*^b2h,^) , . . . ,t.„) h*^. 

Since ti,. . . ,tk are I?- free, the last cumulant vanishes unless ii — ■ ■ ■ — in (and 
hence /iiti/i*, . . . , hkikh^ are S-free). Moreover, in the latter case we have that 

Kn (hitih*b, . . . , hiUh*) = h^K^ {UF {h*bhi) ,.. .,ti)h* 

is C-valued for all 6i, . . . , bk-i G B (in particular for all 6i, . . . , fe^-i G C), and 
thus, by Proposition 12. 1[ the C-cumulants of ft-itift,*, . . . , hkikh^. are the restrictions 
of the S-cumulants. Thus the vanishing of the mixed S-cumulants implies also 
the vanishing of the mixed C-cumulants, hence /iitift,*, . . . , hktkh"^ are also C-free. 
Furthermore, we get, for all c G C, and alH = 1, . . . , fc that 

i^i^i^i ■ • ■ , hitih^ ) — {hitih^ c, . . . , hitih^ ) 

= hiK^ (tiF {h*chi) ,U) h* 

= hiKn {tiPiT^'''\h*chi), ...,U^ h* 
= h,h* (T^'\h*ch,)y K« iU, 

Hence 

oo 

T^hitih^ (c) = {hitih*c, . . . , hitih*) 

oo _^ 

= Y^h^h* {t^^ {h*ch,)Y ^{t„...,t^ 
n=l 

Since hitihl, . . . , hkikh^. are C-free we get for the C-valucd 7^-transform of $ that 

fc 

7^g(c) = ^7^^^,,^J(c). 

i=l 
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Now we apply all this to our functional relation between the C- valued Cauchy and 
7?.-transform. We obtain 



G%{c) = {n%{G%{c))-c) 



1=1 J 
k \ 
Y^h^Knt^ (r« {h*Gi{c)h,))-c 



and therefore 



i=l 



(4.1) = (^E'^^'^^*^^'^ (^^^°^(^r(Pocpo)/ii/i:)) -Poc^^oj , 

where in the second step wc introduced the po in the argument of r*^*^ by using 
Pohi = hi; in the third step we expressed t*^*) in terms of t'"' according to their 
definition. 

Consider now c G C^^\ so that pocpo = c. Then, with the definition 

N 



this becomes 



fj{c) :=^t(°) [Gi°{c)h,h;) 



Gf\c)={J2hih*n^^ {fi{c))-c 



and thus we have for the fj (c) the system of equations 



/.(c) 



N 



.(0) 



^/^./^*7^^f (/.(c)) 



In order to get now the scalar-valued Cauchy transform G^{z) (for 2; e C) of 
our free deterministic equivalent $ = we project the operator-valued Cauchy 
transform to the scalar level: 



where we also put fi{z) := /^(zlgco)). 

If we rewrite all this in terms of the matrices Hi and Tj we arrive finally at 



G$ {z) 



/ k \ ■ 

J2HiH*TZTAfi{z))-zlN] 

\i=l J 
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where the fj{z) are determined by the system of equations {j — 1, . . . , k) 

-1 



h (^) = ^Tr 



r k \ 

^i?,ij;7eT.(/.(z))-z/jv 

\i=l / 



These equations are equivalent to the ones showing up in ^ . (Since there they 
do not use the 7?.-transform of the matrices T,;, this information has to be encoded 
in another set of equations in their approach.) 

One should of course also consider the question whether those equations deter- 
mine the fj{z) uniquely, within a suitably chosen class of functions. This questions 
is answered affirmatively for the present example in [BJ. Let us also point out that 
it is reasonable to expect that such a question (in the form whether the operator- 
valued equation (|4.ip has a unique fixed point in the class of Cauchy transforms) 
can be treated in a suitably general frame by extending the approach from [lOj . 

Remark 4.1. In f6], they consider the seemingly more general problem, with deter- 
ministic rii X rii matrices Ti and random matrices Wi = UiPi which are rectangular 
Ni X Hi, obtained by considering only the first rii columns of a square iV^ x Ni 
Haar-distributed unitary matrix. 

If we let Ti be the inclusion of T in the space of Ni x Ni by just filling up with 
zeros, we have 

U^f.U* = U,P,f,P,U* 
= W^T.W*. 

The convergence in distribution of T w.r.t. ;^Tr is equivalent to the convergence in 
distribution of Ti w.r.t. j^Tr (the moments just get scaled by a constant). Hence 
it will be enough to solve the case when Ni = rii. 

5. Asymptotic P-freeness of Randomly Rotated Matrices 

In this section we will prove Theorem l2.7l For notation, definitions and a detailed 
exposition of the results used in this section, we refer to the last lecture of [T5] . 
where the square case is treated. We will only point out those aspects of the proof 
which differ from the square case. First we will compute the joint 2?-moments 
of Vi,\JT>2V* assuming that the desired asymptotic 2?-freeness holds. Then we 
compute our asymptotic joint moments by the use of the formulas relating the joint 
distributions of the entries of Haar unitary matrices and the Weingarten function. 
We realize that both computations coincide. 

5.1. Computation of Operator-valued moments. Suppose that {A,t) is a 
(pi . . . ,pfc) -rectangular probability space and that ci, . . . , Cp, di, . . . ,dq are simple 
elements such that T>i := (ci, . . . , c,}, T>2 := (di, . . . , dp) are P-free. The goal of 
this section is to express the mixed I?- moments of Vi , T>2 in terms of the individual 
"D-moments of 2?i and 2?2 

By linearity we may only consider monomials. The first observation is that a 
product of simple elements is simple as well. So a general 2?-moment will be of the 
form 

F(i?iCi...A.ao, 

where Dj £ is a product of the diS and Ci is a product of the c^'s. In 

order to have a non-zero moment, we must have that Cj G A'^'^ (i).''0+i)). in the 
following we only consider such products. We denote by Vm '■— {x : r{x) — m} C 
[n]. 

If we consider a totally ordered set 5 = 5*1 U ^2, where Si n 5*2 = and we let 
{Vi, . . . , Vs} = TTi be a partition of 5*1 and {Wi, . . . , Wt} = 7r2 be a partition of 5*2, 



FREE DETERMINISTIC EQUIVALENTS 



15 



we will write tti U7r2 for the partition {Vi, . . . , Vs,Wi . . . , Wt} of S. For a partition 
TT G NC{n) we consider its Kreweras complement, defined as the largest partition 
K{Tr) e NC{T, ...,n) such that tt U if(7r) e 7VC(1,T, ...,n,fi). For a pair of non- 
crossing partitions a < tt G NC{n), we have that a U K{Tr) S -/VC(1, 1, . . . ,n,n). 
We will be interested in pairs a < tt such that 



r'{i) = r(cr(i)),r(i) = r'(ii:(7r)(i - 1)), 



1 < i < n, 



where the arithmetic is performed modulo n and a partition tt is viewed here as the 
permutation in S'„ with cycle decomposition prescribed by the block structure of tt 
(ordered increasingly) . The collection of such pairs of partitions will be denoted by 
NC{n)r.r' C NC{n)'^. For such pairs we define the weight 



\ ( 



\i<a{i) 



I 



n ^(^'-(j+i)) 



ie[«] 

j<K(7r)0) 



Proposition 5.1. Let F(Z?iCi . . .DnCn), r and r' be as above. Then 
¥{D^Ci...D^Cn) 

[Di,...,D„]-T^[Ci,... 



, ,C„] • /i„[cr, tt] 



Proof. We use the operator-valued moment-cumulant formulas 
F(i^iCi...i?„C„) = 



pGAfC(14,...,n,ri) 



By freeness of Vi, 7)2, we see that for a non- vanishing contribution to the sum, 
p must be of the form tt U tt, where tt G NC{n) and tt < K{t:). By fixing one 
TT G NC{n) we may sum over all partitions tt < if (tt), to obtain an expression 
where the sum runs over tt G NC{n) and the sumand consists on evaluating in 
a nested way (indicated by tt U if(7r)) the operator- valued cumulants of the D's 
and the operator-valued moments of the C's. Then, for each fixed tt we use the 
cumulant-momcnt formula to express the cumulants in terms of moments, obtaining 



F(i?iCi...A.C„) 



E 



tUK{tv) 



Ti)[Dl,Ci, . . . ,Dn,Cn] ■ ^fATC K tt] . 



a<TTeNC(n) 



Finally it is easy to observe that we only have a non- vanishing contribution when 
(<T, tt) G NC{n)r.r'- lu thcsc cases we can use formula (|2.ip to get the last expres- 
sion. □ 



5.2. Asymptotic Joint Distribution of Haar Unitary Matrices of Different 
Sizes. We conclude the proof of our theorem by showing that the asymptotic T>- 
moments of 'D^\\] i^'D^^JJ*^^ are computed as in Proposition 15. II 

First we fix (and omit) N for simplicity of notation. One can see that a general 
non-vanishing operator-valued moment is of the form 



F([/,(i)i?(i)c/;,(i)C(i) . . . c/,(„)Z?(")c/;,(„)C(")), 
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where D^J') € ^('■(j)''-'^^')) and C^^) S ^('■'0),'-0+i)) are products of simple elements 
in and T>2^'^ , respectively. Then we have 

_ ^r(l) 1 ^/ (r(l)),(l) (r'(l)) (1) (r(n)) ,(n) (r'(«)) („) 

The only difference from [15] is that the expectation of the product of the entries 
of our unitaries will be factored by the assumption of independence between the 
Haar Unitaries. We can apply Lemma 2.2 from [53] to each of these factors. In the 
language of the permutations considered in [TS] this factorization of the expectation 
will be translated to the conditions on a,/? G Sm that = r^(j), ri = v'^^^y 
i — 1, . . . ,n. We denote by Sn,r,r' C Sn the subset of such permutations. If we 
write am ■— a \v^, Pm '■— /? | v„ , we notice that a'^Pm G -Sy,,, • Then one can show 
that the sum of the right hand side of the last equation is actually 



]lfl JnJ„ 1^12 11 

se[n] m=l 



( \ 

n -S! n 



s6[n] s'G[n] 
Yr(s)=m r'{s')=m 



J2 Tr„[Z?(i),...,7?(")].Tr^-i^[C(i),...,C(")]. n Wg(^™,a„i/?™) 

a,/3eS„^^_^/ m=l 

a,/3eS„_^_^/ m=l 

where 7 G S",! is the long cycle 7 = (1, 2, . . . , n — 1, n), 4t^{a) denotes the number of 
cycles of a and each factor Wg is the Weingarten function on Sv,y^ = S\v„^\- Since 
^#(a~^/3m) = #(Q!~^/3), all further remarks on the order of the contribution 
remain the same as in the square case. 

For a non-vanishing limit, we must have again a, /3 6 SNc{n), a < (3. Such pairs 
of permutations are known to be in one-to-one correspondence with pairs Pa < Pp, 
of non-crossing partitions. One can see that the image of Sn,r,r' under P is exactly 
NC{n)r,r' ■ Hence the sum here can be thought as running over the same objects as 
in Proposition lS.ll For these pairs of permutations, we get the limiting contribution 
of 



m— 1 



rn—1 rn—1 

Where </> is the coefficient of the leading term in the expansion of Wg. It is not hard 
to see that the second product coincides with the weight function uj{Pa, Pi3)r,r' of 
the corresponding partitions. 
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Finally, we observe that if a,/3 € S]yc(n)j ol < P, then a £ S]\!c(n)- One 
shows that 



fe 



m=l m—1 

/ k 



Hence, all the factors appearing in our non-vanishing terms correspond with the 
ones in Proposition 15.11 We conclude that I?^^^ UArI?2^''U^ are asymptotically 
P-free. 
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Appendix A. Extension of Estimates for the Difference of Cauchy 

Transforms 

A.l. Introduction. Roughly speaking, we want to discuss how estimates for the 
difference of Cauchy transforms near infinity (i.e. on a subset := {z E C+; \z\ > 
R} of the upper half-plane C"*" := {z G C; Im(z) > 0} with R > 0) extend to 
estimates near the real line. The approach presented here is motivated by a paper 
of V. Kargin ( |A3j ). In order to understand how powerful the method is, we will 
reformulate his results in a more general framework. The main goal is the following 
theorem: 

Theorem A.l. Let /3,c e (0, 1) and R,T > be given. Then there exist constants 
Too > and 770 > such that for any two probability measures fi and v on the real 
line R with 

sup |G^(z)-G,(z)| <e-'" 
for some m > toq the inequality 



holds, where 



max \G^{z) — Gy{z)\ < exp ( — cm 



Tl \ ■ ^ ^^"^ ( m'^ ) I r Til 

/(m) 2770— -, + [-T,T]. 

1 exp ( - ;;^) 



This will enable us to prove the following result about the qualitative behavior 
of the Kolmogorov distance: 

Theorem A. 2. Let fi be a probability measure with compact support contained in 
an interval [—A, A] such that the cumulative distribution function J-"^ satisfies 

(A.l) \J^f,{x + 1) - Tf,{x)\ < p\t\ forallx,teM. 

for some constant p > 0. Then for all R > and P e (0, 1) we can find > and 
ruo > such that for any probability measure v with compact support contained in 
[—A, A], which satisfies 

sup |G^(z)-G,(z)| Ke-" 

zeA+ 

for some to > toq, the Kolmogorov distance d{p,v) := supj,gjj |7^^(a;) — Jv(x)| 
fulfills 

mP 

A. 2. Proof of Theorem lA.li We use some results of complex analysis on the 
unit disc D := {z G C; \z\ < 1}. For any / G 0(10'), i.e. for any holomorphic 
function / on D with values in C, we define 

M(f,r) := max |/(re''')| for all r G [0, 1). 

eG[0,27r] 

IVToreover, for any given continuous function lo : [0, 1) (0, 00), we consider 
B{u) := {/ G 0(D); Vr G [0,1) : M{f,r)<uj{r)}. 
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The key for the proof of the fonowing theorem is a well-known fact due to G. H. 
Hardy ( [A2j ) . which states that the function 

M(/,-): (-00,0) ^M, ,sK^log(Af(/,e^)) 
is convex for any / G 0(D) with / ^ 0. 

Theorem A. 3. Let uj : [0, 1) — > (0, oo) be a continuous function satisfying 
(A.2) lim s" logfwfe"")) = for all a > 0. 

s->0+ 

Then for all rp G (0, 1), /3 G (0, 1) and c G (0, 1) there exists a constant mo > 
such that any function f G ^(a;) fulfilling the condition M{f,rQ) < e~™ for some 
m > mo satisfies 

M(/,exp(^))<»p(-e„..-»)^ 

Proof. We define the function 

(-cx),0)^M, s log(w(e")) 

and put So :— log(ro). If m is sufficiently large, we have that si := ^ and S2 := ^ 
satisfy the condition sq < S2 < si. Then we consider 

I: M ^ M, u{si) - f 1 - )m. 

si - So ^ si-sq/ 

Since 

EI^^Ll^ 1 and m^fl^'-l^)^^— 1 

si-so l-^ \ si-soJ 

we deduce, using the equations 

1{S2) 1 S2"So/, A «/ S2-S0 



■(isirc(si) 



ml |so|°si-soV' ■ ■ '/ \ si - So 

and (jA.2p with a :— > 0, that 1^1^^ — > — 1 as m — cxa. In particular, this 
implies that we can find mo > such that ^52) < —cm}~^ is fulfilled for all 
m > mo- 

Let m > mo be given and let / G B{uj) be a function satisfying M{f,ro) < e"™, 
where we restrict ourself to the non-trivial case / ^ 0. Because M(/, •) is a convex 
function, it follows from the inequalities 

Z(so)--m>log(A/(/,e^°)) =M(/,so) and ?(si) = 5(si) > M(/, si) 

that also the desired inequality Af(/, S2) < '(52) < ~cm^^^ holds. □ 

Since we are interested in the behavior of Cauchy transforms, which are holo- 
morphic functions living on the upper half-plane, we have to translate the above 
theorem. Therefore we define for every a > a biholomorphic mapping 

I + z 

ipa- © ^ C+, z ia . 

1 — z 

Theorem A. 4. Let lj : [0,1) (0, 00) he a continuous function, which satisfies 
condition (jA.2p . and let aQ > be an arbitrary constant. We consider 

B*{u, ao) := {/ G 0(C+); Va > ao : / o G B{oj)}. 

Let 13, c G (0,1) and R,T > be given. Then we can find constants mo > and 
ryo > such that any function f G S*(w,ao), which satisfies sup^g^+ |/(z)| < e~™ 
for some m > niQ, also fulfills 



max 
te[-T,T] 



1 - exp ( - 45-) 



< exp ( — cm 
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Proof. It is an easy exercise to show that the foUowing equahty holds for all z G 

2a\z\ 



pa{zj - la- -— 



and to deduce that ijja maps the disc D{0,r) with r > to the disc D(ia ^'^^^ , -^zpi)- 
If we put ro :— exp(— i) e (0, 1), it is therefore possible to choose a> gq such that 
ipa{D{0,ro j) C is satisfied. Hence we have M{f o tpa,ro) < sup^g^+ An 
application of Theorem IA.3I shows that, in the case where sup^g^+ 1/(^)1 < e~™ 
holds for some m > toq, we also have 

A/(/oV..,exp(ii^))<exp(-c.n- 

This means that |/| is bounded by exp ( — cm^^^) on the disc ■!/'a(-D(0, r)) with 
r := exp ( — 7^73- )• For an appropriate choice of 770 > (an easy calculation shows 
that each ryo > 2 ° works), the disc D(^iaj^^, contains the whole interval 

irjoj^ + [— T, T] for sufficiently large values of m. This leads finally to the desired 
inequality. □ 

The validity of Theorem lA.ll relies on the fact that we can apply this result to 
functions of the form — with probability measures and on M. 
Since ipa - ^ ^ C"*" gives by extension a homeomorphism fpa\dii>\{i} ■ dD\{l} — > M, 
we deduce that the Cauchy transform of any probability measures /i on M 
transforms in the following way 

G^(^,(z)) = / , ,\ M t) = ^ [ ii^iKL^d/2,(0 for ah z e ID 
jRVa(z)-t 2ai 7aD\{i} z-^ 

where the probability measure Jla is given as the push forward of by (V'a|aD\{i})^^- 
Thereby we used the formula 

=:^ ^^~^^V^^ for aU e e aD\{l} and z G B. 

Since |1 - z| < 2, |1 - ^| < 2 and |z - ^| > 1 - \z\ holds for ah ^ e 90 and z G D, 
we easily see MiGf^ o ^p^, r) < ^ ■ for < r < 1. 

If ji and V are two probability measures on M, we deduce with the above result that 

M{{G^~G,)o^,,r) <-.-^<-^ 

a 1 — r 1— r 

holds for all < r < 1 and a > oq := 4. If we define oj : [0, 1) (0, 00) by 
'^(^) •= 137 J it is an easy exercise to show that u! satisfies condition (jA.2[) . Thus 
we get — Gjy G /B*(cj, oo) and Theorem I A. II follows immediately from Theorem 

A. 3. Proof of Theorem IA.2[ At first, we recall a result of Z. D. Bai (Corollary 
2.3 in |Alj ). which particularly states that there are constants ci,C2 > such that 

d{n,iy)<cii / \Gf,{iri + t)~G^{ir) + t)\dt+- sup \F ^,{x + 1) - T f,{x)\dt 

holds for all probability measures fi and v with supports belonging to the interval 
[—A, A] and for all 77 > 0. If /i satisfies (|A.ip . we apply Theorem lA.il in the case 
T = C2A and for any probability measure v with compact support contained in 
[—A, A], which fulfills sup^g^+ |G^(z) — Gp(z)| < e^™ for some m > toq, we choose 

1 - exp ( - ^) 

1 + exp ( - 
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Hence we get 

('C2A 

\G^{iT]{m) + t) - G^{ir]{m) + t)\dt < 2c2Aexp ( - cm^'^) 



C2A 



and by assumption (lA.ip 

\T^j,{x + t)~ F^{x)\dt < p \t\dt = 16p7j{mf. 

-4jj(m) J —iri(m) 

Since linim^oo m^rj{m) = ^rjo and linim_>.oo Tn^ exp{—cm^^^) =^ 0, it follows that 



d{p, I') < ^2ciC2 Am^ exp ( — cm^ + I6cipm'^ri(rn)j 



holds with an appropriate choice of 9 > (for instance, each Q > Sciprjo works) 
for all sufficiently large values of m. 
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